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Abstract 

Complete descriptions of the Lie symmetries of a class of nonlinear 
reaction-diffusion equations with gradient-dependent diffusivity in one 
and two space dimensions are obtained. A surprisingly rich set of Lie 
symmetry algebras depending on the form of diffusivity and source 
(sink) in the equations is derived. It is established that there exists 
a subclass in 1-D space admitting an infinite-dimensional Lie algebra 
of invariance so that it is linearisable. A special power-law diffusivity 
with a fixed exponent, which leads to wider Lie invariance of the equa¬ 
tions in question in 2-D space, is also derived. However, it is shown 
that the diffusion equation without a source term (which often arises 
in applications and is sometimes called the Perona-Malik equation) 
possesses no rich variety of Lie symmetries depending on the form of 
gradient-dependent diffusivity. The results of the Lie symmetry clas¬ 
sification for the reduction to lower dimensionality, and a search for 
exact solutions of the nonlinear 2-D equation with power-law diffusiv¬ 
ity, also are included. 
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1 Introduction 


Nonlinear diffusion equations of the form 


u t = div (D(u)Vm) 


( 1 ) 


(hereafter u = u(t,x ) is an unknown function, x = (aq ,...x n ), D is an ar¬ 
bitrary nonnegative function of its argument, Vit = (u Xl ,..., u Xn ) and the 
subscripts t, x i ... x n , denote differentiation with respect to these variables) 
have an enormous number of diverse applications and have accordingly been 
widely studied by Lie symmetry method for a long time (see e.g. [TUI for 
n — 1, mm for n — 2,3). 

An alternative class with similar status can be written in the form 


u t = div (D (fl) V«), VL = | Vu| 2 , 


( 2 ) 


and also has a significant number of applications (some of which we note be¬ 


low) but has been much less studied from symmetry point of view (excepting 
the case n — 1). An observation relating to the similar standing of (JT]) and 
([2]) is that the later corresponds to a gradient flow 


du 5L 

dt 5u 


involving the variational derivative of 



while the former corresponds to 


du 1 5L 

dt ^'(u) 5u 


with 



A third equation of similar status is 
w t = '&\ Aw), A = V 2 


( 3 ) 
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which is reduced to (HD, now with D{u) = T" by the non-local substitution 
u = Aw: this relationship is a non-local one, so the two may have different 
Lie symmetry properties. Taking into account this substitution and making 
the calculations, © and (J3J) can be written as the system 

du 5L dw 5L 
dt 8w ’ dt 8u ’ 

where 


L 


T( u)dx 


T(A w)dx. 


Finally, we in this case arrive formally at the equation 


dL 

dt 


, 5L l2 7 _ 

| V —— 1 2 dx. 
ou 


In one dimension equations (JTj) , ([2]) and ([3]) are non-locally equivalent and 
we have the hierarchy (replacing u by v in the second equation) 

u t = {^'{u)) xx , v t = (^'(Uj.)) , w t = ^'(w xx ) 


with u = v x and v = w x . Lie symmetry properties of this hierarchy have 
been studied by mm- However, when source (sink) terms are included, as 
we do below, the corresponding reaction-diffusion equations are inequivalent 
even in one dimension. 

It should be noted that equation (J 2 ]) for n — 2 was proposed by Perona 
and Malik in m as an edge-enhancement model in image processing and 
this represent one of the most prominent applications of (J 2 ]) , stimulating a 
great deal of interest in the image processing community during the last two 
decades. It is commonly proposed that the Perona-Malik equation provides 
a valuable algorithm for image segmentation, noise removal, edge detection, 
and image enhancement (see, e.g., [13]). In ((2]), the unknown function u = 
u(t,x 1 ,^ 2 ) describes the time evolution from an original image, while the 
diffusion coefficient D (hereafter Dq 7 ^ 0) is a nonnegative function of the 
magnitude of local image gradient V« = (u Xl ,u X2 ). In [21], the authors 
proposed two forms for the diffusion coefficient D, namely 

B<n) = exp (-^) ’ ° (si >=0+ <4 > 
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(here D 0 is a positive constant), which are the most common in image mod¬ 
elling. Equation ([2]) with the power-law diffusion coefficient 


D(Q) = Q m , m> 0 


(5) 


has also been used in modelling turbulent gas flows: see [4], for example. 

Rigorous mathematical results for the Perona-Malik equation have been 
obtained mostly under the one-dimensional restriction, i.e. 


u t = (D(u x )u x ) x 


( 6 ) 


and relevant references are given in [22] (see also the earlier papers HUE!)- 
Equation (0 is sometimes called the nonlinear filtration equation and is 
applied in modelling non-Newtonian flows of fluids. 

Here we concentrate ourselves on the investigation of (J2D and its general¬ 
isations by means of the Lie symmetry method. Because their Lie symmetry 
properties essentially depend on the space dimensionality, we examine sepa¬ 
rately the cases n — 1 and n — 2. The Lie symmetries of ([6]) were completely 
described many years ago [2]. However, to the best of our knowledge, there 
are no rigorous mathematical results concerning the Lie symmetry of (J2D in 
higher dimensions. Moreover, its extensions involving a source (sink) nat¬ 
urally arise in applications (see, for example, ra, ra>. Thus, one may 
consider the class of equations 


u t = div (D (Q) Vu ) + Q(u), 


(7) 


where D ^ constant and Q are arbitrary smooth functions of their argu¬ 
ments. Equations of the form (JTJ) also have an obvious gradient-flow inter¬ 
pretation. 

In Section 2, a complete description of Lie symmetries of equation (0 
is performed in the one-dimensional case. In Section 3, the same is done in 
two dimensions, the special case (0 being examined separately. It is worth 
emphasising at the outset that a number of the symmetries obtained are far 
from obvious by inspection. In Section 4, radially-symmetric reductions of 
equations with the power-law diffusivity to ODEs together with examples 
of exact solutions are presented. Finally, we discuss the results obtained 
in order to compare those with the known results derived for more familiar 
nonlinear reaction-diffusion equations. 
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2 Lie symmetry of a class of (l + l)-dimensional 
equations 

In this section, we consider the class of equations 

Ut D(Ux)U xx T Q(w), (8) 

where the smooth function D(u x ) is an analog of the diffusivity in ([ 6 ]) (ob¬ 
viously (ED is equivalent to the (l+l)-dimensional equation (0) on redefining 
D suitably). In what follows we assume that D{u x ) ^ constant because 
(ED with the constant diffusivity is the standard reaction-diffusion equation, 
for which the Lie symmetry classification was performed in [9], It is also 
assumed that Q{u ) ^ constant because a constant source (sink) Q(u ) = q in 
(ED is removable by the substitution 

t = t, x = x, u —>■ u — qt, (9) 

i.e. (ED with Q(u) = q is equivalent to (ED- 

First of all we note that equation (ED with D(u x ) = u~ 2 is linearisable 
by the well-known hodograph transformation. This linearisation is a direct 
consequence of Lie symmetry properties of ((HD with D{u x ) = u~ 2 . In fact, 
by direct application of the Lie algorithm |IZ:,il9] one proves the following 
statement (see a crucial step of the proof presented below after Theorem 3). 

Theorem 1. The nonlinear equation 

u t = u~ 2 u xx + Q{u), (10) 

with an arbitrary smooth function Q{u ) admits an infinite-dimensional alge¬ 
bra of Lie symmetries generated by the operator 

X°° = w(t,u)d x , (11) 

where w(t,u ) is an arbitrary solution of the linear PDE 

wt = w uu - Q{u)w u . ( 12 ) 

The hodograph substitution x —> u, u —* x transforms OB exactly to (fl 2 l) 
with w(t , u ) = x{t , u ). 
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Taking into account the special status of D{u x ) = u~ 2 , we consider fur¬ 
ther equation (jSJ) with D{u x ) 7 ^ u~ 2 . In order to obtain a complete description 
of Lie symmetry operators of equations belonging to class (JHD, we start by 
finding of the group £ of equivalence transformations as the first step. By 
definition, each transformation from £ transforms an arbitrary fixed equa¬ 
tion of the form (J 2 S 1 ) to another equation (which can be the same equation) 
from this class. Nowadays finding the group £ for a class of PDEs is a stan¬ 
dard procedure, which consists in application of a well-known algorithm (see 
e-g- H), hence we present the result only. 

Theorem 2. The group £ of the class of PDEs (JB]) consists of the group £ c 
of continuous equivalence transformations generated by transformations 

~ d 2 ~ 'y 

t — at + 6 0 , x = /3x + 5 1 , u = -yu + <f 3 , D = -— D , Q = — Q, (13) 

a a 

(here a > 0 , (3 > 0 , 7 > 0 , h 0 ; $i and S 3 are arbitrary real constants ) and 
the set of discrete transformations 

1. x —* — x; 

2. t ->■ —t, D ->■ —D, Q ->• —Q; 

3. u —y — u, Q —y — Q. 

Remark 1. The discrete transformation 2 is valid under the assumption that 
time and diffusivity may be negative, which means that the results presented 
below are also valid for the associated non-physical (backward) equations. 

Obviously, the discrete transformations listed above can be incorporated 
into the formulae (TT3T) by allowing a, (3 and 7 also to be negative, i.e. only 
the restriction a(3 7 7 ^ 0 takes place. 

Now we find so called the principal algebra (other terminology used for 
this algebra is the ‘trivial algebra’ and the ‘kernel of maximal invariance 
algebras’) of the nonlinear PDEs class () 8 l) . i.e. the maximal invariance algebra 
(MAI) admitted by each equation of the form (jHJ). A direct application of 
the classical Lie algorithm under the assumption that the coefficients D(u x ) 
and Q(u) are arbitrary smooth functions leads the two-dimensional Abelian 
algebra ALie 2 generated by the basic operators 

X 1 = d t , X 2 = d x . (14) 
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Tabic 1: Lie symmetries of the class of equations (JEJ) with D ^ u x 2 



D(u x ) 

Q{u) 

Constraints 

MAI 

1. 

V 

u~ x 


ALie 2 , X-,i = 2 td t + xd x + ud u 

2. 

V 

u 


ALie 2 , X 3 = e t d u 

3. 

< 

£\e~ u 


A.Lie 2 , = (/c -(- 2 )tdt xd x -K {h 2)<9 W 

4. 

< 

El'll™ 

m / 1,2 

ALie 2 , X 3 = (1 - m)td t + k ff~™ x d x + ud u 

5. 

< 

£\u k+1 + £ 2 u 

kA± i 

ALie 2 , X 3 = e~ ks2t (d t + £ 2 ud u ) 

6. 

1 

U x +J 

u 

7 ¥= o 

ALie 2 , X 3 = e Slt d u , X 4 = e £lt (d t + £1 (u + jx)d u ) 

7. 

U x 

£\U 2 


ALie 2 , X 3 = td t — ud u , X 4 = t 2 d t — (2 tu + £\)d u 

8. 

u x 

e\U 2 + £ 2 

£i£ 2 = -1 

ALie 2 , X 3 = e~ 2t {d t + 2 [u - £i)d u ), 

X 4 = e 2t (d t - 2 (u + £i)d u ) 

9. 

U x 

£\U 2 + £2 

£l£ 2 = 1 

ALie 2 , X 3 = cos(2 t)d t + 2(sin(2t)u + £\ cos(2 t))d u , 

X 4 = sin(2 t)d t — 2(cos(2 t)u — £i sin(2t))d„ 

10. 

< 

£ 2 W 


ALie 2 , X 3 = xd x + (l + f) ud u , 

X 4 = e ~ ke2t ( d t + £ 2 ud u ), X 5 = e E2t d u 


In order to obtain a complete description of Lie symmetries of the PDE 
class (IE]), one needs to find all possible pairs ( D,Q ), for which the relevant 
nonlinear PDE is invariant under a MAI of a higher dimensionality than two. 

Theorem 3. All possible MAIs (up to equivalent transformations) of equa¬ 
tion for any fixed pair of the functions ( D , Q ) are presented in Table 1. 
Any other equation of the form invariant under a non-trivial Lie algebra 
is reduced to one of those given in Table 1 by an equivalence transformation 
from £. 

Remark 2. In Table 1, by rescalling and/or reflecting (see Theorem 2) the 
coefficients £\ = ±1 and e 2 = ±1 without loss of generality. Of course, the 
relevant function Q in real world applications may contain arbitrarily given 
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non-zero coefficients Ai and A 2 instead of these. This remark is valid also for 
Table 2 (see Section 3). 

The proof. Here we apply the Lie-Ovsiannikov approach (the name of 
Ovsiannikov arises because he has published a remarkable paper in this direc¬ 
tion, [18]) of the Lie symmetry classification, which is based on the classical 
Lie scheme and a set of equivalence transformations of the differential equa¬ 
tion in question. It should be noted that there are some recent approaches, 
which are important for obtaining the so called canonical list of inequivalent 
equations admitting non-trivial Lie symmetry algebras and allow the solution 
of the problem of the Lie symmetry classification in a more efficient way then 
a formal application of the Lie-Ovsiannikov approach (see some remarks at 
the end of this section). 

Thus, in order to fold Lie symmetry operators, one needs to consider each 
equation of the form ([ 8 ]) (clearly this should be viewed a class of PDEs be¬ 
cause the functions D and Q are arbitrary smooth functions) as the manifold 


M = {u t - D(u x )u xx - Q{u) = 0} (15) 

in the prolonged space of the variables: t, x, u, u t , u x , u xx , u xt , u tt . According 
to the definition, (JHJ) is invariant under the transformations generated by the 
infinitesimal operator 


X = £°(£, x, u)d t + ^(t, x , u)d x + rj(t, x, u)d u , 
if the following invariance condition is satisfied: 


X 


(ut - D(u x )u xx - Q(u )) 


= 0. 


M 


(16) 


(17) 


The operator AT i s the second prolongation of the operator X, i.e. 


X 

2 


X 


Pt 


_d_ 

du t 


_d_ 

du. 


+ +<7z 


d 


du,. 


+ CTtt 


_d_ 

du t 


+ &tx 


d 


du. 


tx 


where the coefficients p and a with relevant subscripts are expressed via the 
functions and 77 by well-known formulae (see, e.g., H3H9]). 

Now we present the so called system of determining equations (DEs), 
obtained by direct calculations using the invariance condition (TTT1) 


(Vx + 0 lu ~ Qu x - 6X) D ' + (ft - 2 £x - 2 6 x) D = 0, (18) 
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~Vt + (£°— Vu)Q + vQ + (£* + £lQ) u x+ 

(vxx + (2 v xu - W + {vuu - 2£lu) u x - ZL u l) D = ° ( 19 ) 

(hereafter D' = J-j-p and Q — together with the standard restriction 
£° = £°(i), which is common for any evolution second-order PDE. 

In contrast to the system of DEs for the standard RD equation (see those 
in [9]), which immediately specifies the structure of the unknown functions 
^°,^ 1 and ? 7 , DEs R 8 ll and R91) are much more complicated because the 
diffusivity D depends on u x (other than on u ). The crucial step is to solve 
(1I5D w.r.t. D, taking into account that the functions and rj do not 

depend on u x , so that one may formally write (|T 8 j) as follows 

(eo + u x — e2 u 2 x )D' + [e^ — 2e2 u x )D = 0, (20) 

where e^, k — 0 ,..., 3 correspond to the relevant coefficients in (fT 8 ]l . 

It has been noted that there exists a unique case D = Cu ~ 2 (the constant 
C can be removed by scaling, see Theorem 2 ) when the parameter e 2 can 
be arbitrary in (l20l) because only in this case does the relationship e 2 u 2 x D' = 
2 e 2 u x D apply. Thus, (1201) with D = Cu~ 2 is equivalent to eo = 0, = 2ei, 

i.e. 


Vx = 0, 4° = 2r] u 

while is still an arbitrary function. The remaining DE (TT9|) with D = Cu~ 2 
and an arbitrary smooth function Q has the solution = w(t,u), where 
w(t,u ) is the general solution of the linear equation (TT 2 |) . This means that 
MAI of (Roll is infinite-dimensional (see operator RT]) ) As a result, Roll can 
be linearised (see Theorem 1). Thus, we assume D ^ Cu~ 2 in what follows. 

Solving the linear ODE (l20ll . one arrives on five essentially different cases 
depending on the parameters. 

(i) If e 2 = e\ = eo = 0 then immediately e 3 = 0 and D is an arbitrary 
non-constant function. 

(n) If e 2 = 0 and e\ ^ 0 then the general solutions is 

e 3 

D = c(u x + ^-y^, e 3 ^ 0 . ( 21 ) 
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(Hi) If e 2 = e\ = 0 and e 0 e 3 7 ^ 0 then the general solution is 


D = Ce E o’ 


(iv) If e 2 7 ^ 0 and e\ + 7 ^ 0 then the general solutions is 

f 2 e 2 M x e 3 


D = C exp 


9 du x 

eiu% - ei-u* - e 0 


(v) If e 2 7 ^ 0, ei = e 0 = 0 and e 3 7 ^ 0 then the general solutions is 


D = Cu~ 2 exp ( — —u~ l ). 

Consider case (%). Having e 3 = e 2 = e\ = e 0 = 0 and substituting instead 
of e*;, k — 0,..., 3, the relevant coefficients from (1181) . one arrives at a linear 
system of PDEs with the general solution 

f‘ = ^P + Mi), = + ( 22 ) 

where £°(t),b(t) and p(t) are arbitrary functions. Substituting (122|) into DE 
m, one immediately derives a linear ODE w.r.t. Q(u) because the sec¬ 
ond line of (fTUj) simply vanishes. If Q(u) is an arbitrary function then only 
the principal algebra (fl4l) occurs. Non-trivial results are obtained only for 
Q(u) = An -1 and Q(u) = Xu. In both cases, the arbitrary non-zero constant 
A can be removed from the relevant equations using the equivalence trans¬ 
formations from Theorem 2. Thus, the cases 1 and 2 in Table 1 have been 
obtained. 

Consider case (n), which is the most complicated. Because the function 
D depend only on u x , one derives from (l2Tj) the equalities e 3 = —ke 1 7 ^ 
0 , eo = yei and e 2 = 0 (hereafter k 7 ^ 0 and 7 are arbitrary constants). 
Taking into account the relevant coefficients from (|18p . one obtains the linear 
system of PDEs 


= -KVu-Q, Vx = ^(Vu~Q, £i = o. (23) 


There are two essentially different subcases (iia) 7 7 ^ 0 and (iib) 7 = 0 . 

In subcase (iia), equations (123|) immediately produce = 0 provided 
k 7 ^ —2 (this special exponent was checked separately; nothing interesting 
was found, however). Having £ xx = 0 and equations (I25j) . one notes that DE 
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(fTUl) again reduces to a linear ODE w.r.t. Q(u). Moreover, only the functions 
Q(u) = Ait -1 and Q(u) = Xu lead to non-trivial Lie symmetry algebras while 
the diffusivity is D — C(u x + y) fc . However, these algebras are nothing else 
but those arising in case 1 and 2 of Table 1 provided k ^ — 1. The special 
exponent k = — 1 was examined and the case 6 of Table 1 was found (the 
constants C and A are removable using the equivalence transformations from 
Theorem 2). 

In subcase (iib) 7 = 0, equations (123]) take form 
ft° = (2 + k)£ - krj u , rj x = 0, £ = 0, (24) 


so that = 0 ( k 7^ — 2 otherwise D = Cu~ 2 ). Thus, the remaining DE 
(jT9|) again reduces to a linear ODE w.r.t. Q(u). Its solution for arbitrary 
k t^— 2 leads to the cases 3-5 and 10 of Table 1. Moreover, it was found that 
case 5 is not valid provided k = ± 1 (see this restriction in Table 1), because 
that for k — — 1 is simply a subcase of case 4, while that for k = 1 leads to 
the very unusual Lie symmetry algebras presented in the cases 7-9 of Table 
1 . 

Finally, it was proved by the examination of cases (iii)-(v) that there 
are not any further equations of the form () 8 j), which are invariant w.r.t. a 
MAI of dimensionality three or higher. 

The proof is now completed. ■ 

Thus, we have proved that there are 11 distinct cases in which equation (j 8 j) 
admits a non-trivial MAI, i.e. is invariant under three- or higher-dimensional 
Lie algebra. The case of the diffusivity D(u x ) = u~ 2 is exceptional because 
the MAI is infinite-dimensional, hence, the corresponding equation is linearis- 
ablc independently on the form of the function Q(u). The other 10 cases are 
presented in Table 1. It is noteworthy that this classification is surprisingly 
rich: several of these Lie symmetries (see, e.g. cases 6-9) would be unlikely 
to be identified by using the known results for the standard reaction-diffusion 
equation [9]. 

Finally in this context, we present the following observation. The substi¬ 
tution 


r = — -e £2kt , x = x, w — e £2t u 
£2 k 


(25) 


reduces the equations listed in cases 5 and 10 to the equation with m — k + 1 
listed in case 4 and to the equation w T = w k w xx , respectively. Thus, Table 1 
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can be shortened to 8 cases. However, substitution (I25|) does not belong to the 
group £, hence we have found so called form-preserving transformations [IB] 

(it should be noted that this kind of transformations was introduced earlier 
in [15] but without above terminology) allowing a reduction of the number 
of cases obtained via the classical Lie-Ovsiannikov algorithm (see extensive 
discussions in 0. ra. 0 on this matter). 

3 Lie symmetry of a class of (l+2)-dimensional 
equations 

Here we consider a class of equations of the form 

ut = div(D(ft)Vu) + Q(n), ft = |Vw| 2 , (26) 

where u = u(t, Xi, x 2 ) G M, t G M, (xi, x 2 ) G M 2 ; £>(ft) and Q{u) are arbitrary 
smooth functions of their arguments. It can be easily noted that equation 
(|26|) with Q{u) — q = constant again reduces to (J 2 ]) , by the substitution 

t — t, x\ — x±, x 2 = x 2 , u — y u — qt, (27) 

so that we assume in what follows that either Q(u) is non-constant or Q(u) = 

0 . 


3.1 The case of a non-constant source 

We start from the general case when Q(u) is an arbitrary non-constant 
function. Similarly to the (l+l)-dimensional case, we find the group £ of 
equivalence transformations as the first step. 


Theorem 4. The group £ c of continuous equivalence transformations of the 
class of PDEs ([Mil consists of the transformations 


t = at + 6 0 , x a = /3-f ab x b + 6 a , u = y-u + S 3 , D 



a 


~Q, (28) 

a 


where a, j3, 7 , 5^ (k — 0,..., 3), 7 ab (a, b = 1,2) are arbitrary real constants 
satisfying the conditions a > 0, fd > 0, 7 > 0, and the matrix (7 ab ) G 50(2). 


It can be noted that the set of discrete transformations 
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1 . x\ —y —x\, x 2 —> — x 2 ; 

2. t ->■ —t, D -> - D , Q ->■ -Q; 

3. u —y —u, Q —y — Q. 

also belongs to the group £, hence, one may claim that the group of equiva¬ 
lence transformations of the nonlinear PDEs class (l26j) consists of the trans¬ 
formations defined by formulas (1281) with arbitrary real non-zero parameters 

a, /3, and 7 . 

As the second step, we find the principal algebra i.e. the MAI admitted by 
each equation of the form ( 1261 ) . which comprises the operators of translations 
and rotations. 

Theorem 5. The principal algebra of the PDEs class (128|) is the four¬ 
dimensional Lie algebra ALie 4 generated by the basic operators 

X\ <9^, X 2 i X 3 <9x2) X 4 ^2<9xi xid X2 ■ 

The proof of this theorem consists in the direct application of the classical 
Lie algorithm under the assumption that the coefficients D(£l) and Q{u ) are 
arbitrary smooth functions. The relevant calculations are omitted here be¬ 
cause they are simple but cumbersome and the result (namely the operators 
of translations and rotations) is to be expected. 

The third and final step is highly nontrivial: to fold all pairs ( D,Q ), 
when the relevant nonlinear PDE (1261) is invariant under a MAI of a higher 
dimensionality than four. Note that we present the result obtained in two 
tables because the diffusion coefficient D = fl ” 1 is special and this case, 
which corresponds to a limit case of the second formula in 03]), is especially 
noteworthy. 

Theorem 6. All possible MAIs (up to the equivalence transformations) of 
equation for any fixed pair of the function ( D , Q ) are presented in Tables 
2 and 3, where the designations £1 = ±1, e 2 = ±1 are used and k 7 ^ 0, A 7 ^ 0, 
m are arbitrary real constants. Any other equation of the form LEW invariant 
under a non-trivial Lie algebra is reduced to one of those given in Tables 2 
and 3 by an equivalence transformation from the group £. 

Remark 3. Similarly to the (l+l)-dimensional case, the substitution 

r = t 77 —re 2£2fct , aq = x u x 2 = x 2 , w = e~ £2t u (29) 
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Table 2: Lie symmetries of the class of equations (126|) : 1 



D(Q) 

Q(u) 

Constraints 

MAI 

1. 

V 

u~ x 


ALie 4 , X 5 = 2td t + x 4 d Xl + x 2 d X2 + ud u 

2 . 

V 

U 


ALie 4 , X$ = e t d u 

3. 

n k 

£\e~ u 


ALie 4 , X 5 = 2 (k + 1 )td t + xid Xl + X 2 d X2 + 

2 (k + l)d u 

4. 

n k 

£l U m 

to 7 ^ 0 , 1 , 2 

ALie^, X 5 — 2 tdt + m - 2 k-i+k(m+i) x 





(*1^1 + *2§r a ) m—2/c—l+fc(m+l) 

5. 

n k 

£\U 2k+1 + £2 U 

b + _1 I 

ry 2 ’ 2 

ALie 4 , X 5 = e~ 2 ke 2 t (d t + £2 ud u ) 

6. 


£lU 2 


ALie 4 , X 5 = td t — ud u , A 6 = t 2 d t — (2 tu + £\)d u 

7. 


£l U 2 + £2 

£l£ 2 = -1 

ALie 4 , X 5 = e~ 2 t (d t + 2(u - £i)d u ), 

X e = e 2 t (d t - 2(u + £i)d u ) 

8 . 


£l U 2 + £ 2 

£l£ 2 = 1 

ALie 4 , X 5 = cos(2t)9 t + 2(sin(2i)u + £1 cos(2 t))d u , 

X 6 = sin(2 t)dt — 2(cos(2t)zi — £1 sin(2t))d u 

9. 

n k 

£2 U 


ALie 4 , X 5 = xid Xl + x 2 d X2 + (l + £) ud u , 

A 6 = e~ 2kE2t (d t + £ 2 ud u ), X? = e £ 2 t d u 


Table 3: Lie symmetries of the class of equations (126)) : D = Q 1 



Q(u) 

MAI 

1 . 

V 

ALie 4 , AT = xid Xl + x 2 d X2 

2 . 

A u ^ 1 + £2 u 

ALie 4 , X 5 = x\d Xl + x 2 d X2 , X 6 = e 2 e 2 t (d t + £ 2 ud u ) 

3. 

A u -1 

ALie 4 , -^5 — ^ 1 ^x 1 •£ 2 ^x 21 -^6 — H - udu 

4. 

Xu 

ALie 4 , A 5 = a"i<9 Xl + a; 2 9x 2 

A 6 = e 2 Xt (d t + Xud u ), X 7 = e xt d u 
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reduces the equations listed in cases 5 and 9 to the equation with m = 2k +1 
listed in case 4 and to the equation w T = div (| V'in| 2fc Vw), respectively. 
Remark 4- Substitution (t29|) with k = —1 reduces the equations listed in 
cases 2 and 4 of Table 3 to the equation listed in case 3 and to the equation 
w T = div (| Vw| _ 2 Vr(;), respectively. Thus, Table 3 consists of two essentially 
different cases only; however, both cases have no analogs in Table 2. 


3.2 Nonlinear diffusion equation with gradient-dependent 
diffusivity 

It is quite natural to present the result for the nonlinear equation (J2D 
(hereafter with n = 2 ) separately because of its importance for direct appli¬ 
cations. 


Theorem 7. The group £ of equivalence transformations of the class of 
equations ([ 2 ]) consists of the group of continuous equivalence transformations 


t = at + 5 0 , x a = (3-f ab x b + 6 a , u = ju + 5 3 , D 



a 


(30) 


(the restrictions on the parameters are the same as in (|28j) ) and the discrete 
transformations 


1.x i —>■ —X\; 


2. t —» — t, D —> —D; 

3. u —» — u. 


Similarly to the general class (J2H]), the group of equivalence transforma¬ 
tions of (J2]) is represented only by the formulae (1301) . though with arbitrary 
non-zero parameters a, /3, and 7 . 

I 11 comparison to the principal algebra of the PDE class (l26]h that of ([2]) 
is more complicated but still can be easily calculated by direct application 
of the Lie algorithm, hence the following statement has been proved. 

Theorem 8. The principal algebra of ([2]) is the six-dimensional Lie algebra 
ALiee generated by the basic operators 

A-2: X 3 , X 4 . X 5 = d u: X 6 = 2td t + X\d xi + X 2 d X2 + ud u (31) 

(see X\. X 2 , X 3 , and X 4 in Theorem \Ef. 


15 


Remark 5. The Lie algebra ALie 6 is a semidirect sum of the extended 
Euclidean algebra with the basic operators X 2 , X :i . X 4 , X 6 and the Abelian 
algebra with the basic operators Xj , X 5 . 

As the Enal step, we perform an analysis in order to End all possible 
functions D(Q) leading to extensions of the Lie algebra ALie It turns out 
that the result is (perhaps surprisingly) rather simple. 

Theorem 9. Equation (J2]) admits a MAI of a higher dimensionality than six 
if and only if the function D(Q) = Q k = |Vu| 2fc ; k 0 (up to the equivalence 
transformations from the group £). The relevant MAI is seven-dimensional 
and generated by the operators and the operator 

X 7 = 2 ktd t — ud u . (32) 

Corollary 1. Equation ([2]) with the diffusion coefficients D of the form ((4j) 
is invariant with respect to the MAI with the basic operators (T5T1) . 


4 Reductions of the nonlinear equation (M) to 
ODEs 


According to Theorem 9, equation ([2]) possesses its widest Lie symmetry 
in the case of the power-law diffusivity. Moreover, this case is important 
from the potential applicability point of view (see, e.g. ED- Thus, here we 
examine the nonlinear (l+2)-dimensional equation 

u t = div(|Vu| 2fc Vu) . (33) 

Using the seven-dimensional Lie algebra with the basic operators (l3T]b 
(l32i) a wide range reductions of (133|) to (l+l)-dimensional PDEs and to 
ODEs can be obtained. Here we restrict ourselves to an important case 
when (fool) describes a process with radial symmetry w.r.t. to space variables 
(compare for instance an example in |20j). This allows to reduce equation 
(f33|) to the form 

U t = \ «‘ +1 ) r (34) 

(here we assume U r > 0, while the case U r < 0 leads to U t = — \ (r(— U r ) 2k+1 ) r , 
which can be treated in the same way) by the ansatz 


u(t, X\, x 2 ) 


U(t,r), 



(35) 
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generated by the operator of rotations A" 4 = x 2 d xi — x±d X2 . 

First of all we note that equation (I34p with k = —1, i.e. 

Vt = ; 0 rU D, (36) 

is transformed to the linear equation 

V t = -V zz (37) 

by the hodograph type substitution 

r = W, U — z. (38) 

Thus, we assume in what follows that k ^ —1. Moreover, we assume that 
k 0 (otherwise the linear radially-symmetric heat equation is obtained) 
and k ^ —\ because equation (1541) degenerates to a linear ODE in this case. 
Remark 6. In order to describe a diffusion process equation (j34f) should be 
taken with the sign of 2/c+l on the RHS. Especially taking k — — 1 one arrives 
at the linear heat equation instead of (1571) . However, this observation does 
not affect any results presented below because both equations are equivalent 
up to the discrete transformation t —> —t. 

It can be easily calculated that the MAI of (134)) is four-dimensional Lie 
algebra with the basic operators 

Xo = du, X 1 = d t , D 0 = 2(k + l)td t + rd r , Di = krd r + (k + l)Udu (39) 

i.e. the reduction of (1551) to the radial form partly inherits the relevant 
symmetries of this equation (see (1511) and ([52]) ) and no others. A simple 
analysis shows that this Lie algebra leads to four essential different reductions 
of (1551) to ODEs (any other reduction is a composition of one of these four and 
the relevant group of transformations generated by Lie algebra (159]) ). Four 
essentially different combinations of the operators from (1591) . the relevant 
ansatze and ODEs are as follows. 

Case (i) X \, U = 4>{r). ODE for 0(r) : 

K fc+1 ) r = 0. (40) 

Case (ii) D 0 + AA 0 , U = 2 (k+i) ^ + 0( a; )) ^ = rt 2 < fc+1 ). ODE for 

: 

(2* + w' ++ WZT)^ 1 -* - =0- (41) 
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Case (m) Di + AX], U = e * } 0(a;), ui = re ^, A 7 ^ 0. ODE for 4>(co) : 

(2k + 1)0” + u+V + Tu((j)) 1 ~ 2k - ^ \ 1 0(0')~ 2A ' = 0. (42) 

A A 

fc+i 

The subcase with A = 0 leads to the separable ansatz U = r~k~ 0(f), which 
produces the first-order ODE 


3fc + l\ /k + l\ 1 + 2fe 


k 


k 


a+2k 


(43) 


Case (iv) D 0 + A Di, A ^ 0, U — u> = rt 7 , 7 

0M : 


1+Afc 
2(fc+l) * 


ODE for 


(2k + 1 ) 0 ” +oj V + 7 u;( 0 / ) 1 2k - ~ 0 ( 0 / ) 2A = 0 . 


(44) 


It is a simple task to integrate the nonlinear ODEs (|40l) and ()43l) because 
of their simplicity and we omit this. The nonlinear ODEs (TUT) . (1421) and (1441) 
are not integrable for arbitrary parameters (i.e. similarity exponents). 

ODE (1411) is equivalent to the first-order ODE 


(2k + l)y + u y + 


. 1—2 k 


A 


2 (k + 1 ) 


2(k + 1 ) 


y- 2k = 0 , 


(45) 


where y(cu) —(/)'■ To the best of our knowledge the general solution of 
ODE (144]) for arbitrary parameters k and A is unknown (at least in terms of 
elementary functions). However, we note that this equation for A = 0 is a 
Bernoulli ODE and can hence be linearised. Hence the general solution of 
ODE (144]) with A = 0 of the form 


r / —koj 2 + (7ia; -2fc /( 2fc+1 ) \ h 

J l 2(k + l)(3k + 1 ) ) 


dco + C 2 


(46) 


is constructed (hereafter C\ and C 2 are arbitrary constants). Thus, one 
obtains a two-parameter family of exact solutions for the nonlinear ( 1 + 2 )- 
dimensional equation (144]) of the form 


u = 


f /—koj 2 + C\UJ~ 2k ^ 2k+l \^k 

J V 2(k + l)(3k + 1 ) ) 


du + C 2 , 


OJ = t 2 CfcTTy 
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Obviously, the integral in the RHS of (PX71) can be in terms of elementary 
functions in some cases, for example, for C\ — 0 and for any C i ^ 0 provided 

u ~ I I I 

2 ’ 4 ’ 6 ' ' ' ' 

We next consider ODE (1421) . which, unlike to the previous case, cannot 
be reduced to a first-order ODE for arbitrary parameters k and A. However, 
this equation can be rewritten in the form 



and one may note the special case k + 1 = —2k, i.e. k — — |, when this 
equation is equivalent to the first-order ODE 

uj(cj )) s = — u) 2 (j) + Ci. 


This ODE with C i ^ 0 has no solutions in terms of elementary functions 
while setting C\ — 0 one easily finds 


0 =± 


2 

yc 2 - 2(3A)-C S ' 


(48) 


Thus, one obtains exact solutions of the nonlinear (l+2)-dimensional equa¬ 
tion (1331) with k = — A of the form 


u — ± 



2 

2(3A)- 3 (xf + xl) 2 


(49) 


Finally we turn to the nonlinear ODE (l44|h A detailed analysis of this 
equation was undertaken in [4] (see equation (3.6) with k — 1 therein). In 
particular, it was noted that there is special case 7 = — |, i.e., A = — 3^7 
(see (Hg) above) when this equation can be reduced to the first-order ODE 


u{4>) 


2k+l 


1 


2(3fc + 1) 


+ Ci 


in the same way as we have done above for ODE (1421) . The exact solution 
obtained (see formulae (2.3) in [4]) was used for solving the Cauchy problem 
with the Dirac function as the initial profile. 
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5 Discussion 


In this paper, complete descriptions of Lie symmetries of equation (J7J) 
in 1-D and 2-D spaces (i.e. one and two space variables, respectively) have 
been obtained. As one may note, the results depend essentially on the space 
dimensionality of ((TJ) (so that 3-D and higher-dimensionality cases should be 
also examined) and on the form of the pair (D,Q). We have established a 
special case (see Theorem 1) when this equation in 1-D space possesses an 
infinite-dimensional Lie algebra of invariance so that it is linearisable. We 
also have found the special diffusivity D = |Vu|~ 2 , which is special one of (J7J) 
in 2-D space (see Table 3). It should be stressed that equation (J7J) depending 
on the form of (D, Q) admits a surprisingly rich set of Lie symmetry algebras 
(see Tables 1 and 2). Several of these Lie symmetries would be unlikely to 
be identified by a simple extrapolation from the known results 0HO] for the 
standard reaction-diffusion equation 

u t = div (D(u) V«) + Q{u). 

We have examined separately equation ([2]) with n = 2, which sometimes 
is called the Perona-Malik equation, because of its importance for direct ap¬ 
plications. In contrast to the (l+2)-dimensional nonlinear diffusion equation 

(HD, i.e. 

u t = div (D(u)'Vu ), 

which admits a wide range of MAIs depending on the form the diffusivity 
D(u) |l2j . this alternative equation admits a unique extension of the MAI 
(see Theorem 9). In this sense equation ((2D is not analogous to the standard 
nonlinear diffusion equation (HD (at least from the Lie symmetry point of 
view). Moreover, equation ([2D with the common nonlinearities ([4D does not 
lead to any extension of the Lie symmetry. We believe it plausible that a 
search for conditional symmetries of the nonlinear (1+2)-dimensional equa¬ 
tion (J2D in order to compare with the known result for the (l+2)-dimensional 
nonlinear diffusion equation © 0 could lead to more optimistic results. 

On the other hand, one may consider the (l+2)-dimensional equation 

u t = D(Q)'V 2 u , (50) 

instead of (J2D • I 11 contrast to the (l+l)-dimensional case, equation (l50P is 
not equivalent to (J2D- On the other hand, this equation possesses an infinite¬ 
dimensional Lie algebra in the case D(fl) = D _1 . Moreover this algebra 
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contains the operators of the form 

X°° = A(x i,x 2 )d xi + B(x 1 ,x 2 )d xa 

where the pair (A(xi, x 2 ), B(xi, x 2 )) is an arbitrary solution of the Cauchy- 
Riemann system, having similar structure to the operators 

X°° = A(x 1 ,x 2 )d xi + B(x 1 ,x 2 )d X2 - 2 A xi ud u 

that arise in the MAI for the nonlinear diffusion equation ([!]) with D(u) = 
u~ l . Thus, investigation of Lie symmetry properties of equation (l50j) in 
higher dimensions may lead to essentially different results comparing with 
©• 

Finally, we have reduced the (l+2)-dimensional equation ([2]) with a power- 
law diffusivity, i.e. equation (155]) . to a corresponding (l+l)-dimensional 
equation in order to hnd exact solutions. It turns out that the equation 
obtained still possesses a rich symmetry, hence further reductions to ODEs 
have been performed. As a result, two families of new exact solutions of equa¬ 
tion (1521) have been constructed (see formulae (1571) and (1551) ) and solutions 
found much earlier in [5] have been recovered. 
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